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We present measurements of the semi-inclusive cross-sections for and i/p-nucleon deep in- 
elastic scattering interactions with two oppositely charged muons in the final state. These events 
dominantly arise from production of a charm quark during the scattering process. The measurement 
was obtained from the analysis of 5102 u^-mA\iceA and 1458 t'^j-induced events collected with the 
NuTeV detector exposed to a sign-selected beam at the Fermilab Tevatron. We also extract a cross- 
section measurement from a re-analysis of 5030 v^-'mAnceA and 1060 i/^-induced events collected 
from the exposure of the same detector to a quad-triplet beam by the CCFR experiment. The re- 
sults are combined to obtain the most statistically precise measurement of neutrino-induced dimuon 
production cross-sections to date. These measurements should be of broad use to phenomenologists 
interested in the dynamics of charm production, the strangeness content of the nucleon, and the 
CKM matrix element Vcd- 

I. INTRODUCTION 



Oppositely-charged dimuon production in neutrino-nucleon deep inelastic scattering (DIS) provide an excellent 
source of information on the structure of the nucleon, the dynamics of heavy quark production, and the values of 
several fundamental parameters of the Standard Model (SM) of particle physics. These events are produced most 
commonly in charged-current (CC) neutrino DIS interactions when the incoming neutrino scatters off a strange or 
down quark to produce a charm quark in the final state, which subsequently fragments into a charmed hadron that 
decays semi-muonically. This distinct signature is easy to identify and measure in massive detectors, which allows for 
the collection of high statistics data samples. Consequently, dimuon events have played a significant role in the last 
20 years in understanding charm production in DIS |l|]. 

Since the contribution of the down quark to charm production is Cabibbo suppressed, scattering off a strange 
quark is responsible for a significant fraction of the total dimuon rate: ~ 90% in u^N and ~ 50% in v^N scattering. 
Therefore the measurement of this process provides an excellent source for the determination of the strange quark 
parton distribution function (PDF) of the nucleon. Other methods involving measurements of the difference in parity 
violating structure functions from v and v scattering off an isoscalar target and the ratio of parity conserving structure 
functions from charged lepton and neutrino DIS suffer from theoretical and experimental systematic uncertainties. 
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Its prominent role in DIS charm production makes the s quark PDF an important ingredient in tests of two-scale 
quantum chromodynamics (QCD) [^-Q, the two scales being the squared invariant momentum transfer and the 
charm mass nic- The s quark PDF also enters into background calculations for new physics searches at the Tevatron 
and the LHC. For example, in some scenarios the super-symmetric top quark t is best searched for via its loop decay 
to a charm quark and a neutralino 0, t — > cx'^. The dominant background to this mode is from the gluon-strange 
quark process gs —^ W'^c, with the W~ decaying leptonically. Accurate separate measurement of s and s PDFs may 
also shed light on the interplay between perturbative and non-perturbative QCD effects in the nucleon ||^. 

In the following we present a new measurement of dimuon production from scattering of and Dp, beams from an 
iron target in the NuTeV experiment at the Fermilab Tevatron. This measurement exploits the high purity and 
beams produced by the Fermilab Sign-Selected Quadrupole Train (SSQT) to extend acceptance for muons down to 5 
GeV/c momentum without any ambiguity in the determination of the muon from the primary interaction. Our analysis 
proceeds as follows: We first perform a fit to NuTeV and D^^ dimuon data using a leading order (LO) QCD cross- 
section model to obtain values for effective charm mass rric, parameters describing the size and shape of the nucleon 
strange and anti-strange sea, and parameters describing the fragmentation of charm quarks into hadrons followed by 
their subsequent semi-muonic decay. The LO parameters allow us to make contact with previous measurements and 
provide an accurate description of the dimuon data. We then go on to use the LO model for acceptance and resolution 
corrections in extracting, for the first time, the cross-sections da {v^/v^Fe — s- ^^fi^X) /dxdy for forward secondary 
muons tabulated in bins of neutrino energy E,j, Bjorken scaling variable x, and inelasticity y. These cross-section 
tables provide the most model-independent convenient representation of i'^ and dimuon data. They may be used 
to test any cross-section calculation for dimuon production from iron, provided that the model is augmented with a 
fragmentation/decay package such as that provided by PYTHIA We then repeat the cross-section extraction on 
an older CCFR data set |jl^ from the same detector; and, after demonstrating consistency between the CCFR and 
NuTeV data, we combine the cross-section tables to obtain the most statistically precise high energy neutrino dimuon 
production collected to date. We finish by performing a LO QCD fit to the combined cross-section tables and extract 
precise determinations of mc, the CKM parameter Vcd, and size and shape parameters for the strange sea. 



II. LEADING ORDER CHARM PRODUCTION 



Dimuon production from charm depends on three different components: the charm production cross-section, the 
fragmentation of the charm quark to a charmed hadron, and the semileptonic decay of the charmed hadron. In LO 
QCD, charm production arises from scattering off a strange or down quark: + N — > fj,~ + c + X, where 
the second muon is produced from the semileptonic decay of the charmed hadron (Fig. |^). In this case, the dimuon 
production cross-section factorizes into the form: 

d3a^O(.^.jV^^-^+X) ^ d^a^O{.,N^cX) ^^^^ 
dx dy dz rf^ dy 

Here, y is the inelasticity, z is the fraction of charm quark momentum carried by the charm quark hadron, and 
^ = a; (l -I- rn^/Q^^ (l — a;^M^/(3^) is the fraction of the nucleon's momentum carried by the struck quark, with x 
the Bjorken scaling variable and M the nucleon mass. D{z) is the fragmentation function for the charm quark; Be, 
the semi-muonic branching fraction for charmed hadrons; both averaged over all charmed hadrons produced in the 
final state. The charm production cross-section da{vi^i_N cX)/d^dy for an isoscalar target N can be expressed as 

d^a^° {upN fi-cX) _ 2GIME/t: 



didy {l + Q'^/Ml.fV 2ME,( 



1 (2) 

W) ^ 

where Gp = 1.17 x 10"^ GcV^, Mw = 80.42 GeV/c^, and Vcs/Vcd are the cs/cd CKM matrix elements. The 
corresponding i?^ process has the quarks replaced by their anti-quark partners. One observes that the non-strange 
contributions to the cross-section are large in neutrino mode, where they dominate at high ^; and small in anti-neutrino 
mode, where s Q^) dominates at all ^. In the case where s Q^) = s Q^), the neutrino ^ distribution deter- 
mines the relative size of the s Q^) and (u Q^) -f d Q^)) /2 contributions, the anti-neutrino ^ distributions 
determines the shape of s the energy dependence of the cross-section determines mc, the ratio of the dimuon 

cross-section to the single muon inclusive cross-section sets i?c, and the energy distribution of the charm decay muon 
constrains D (z). 
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The non-strange PDFs are determined from inclusive vN and vN CC scattering, and from charged lepton scattering; 
they contribute negligible uncertainty to <Tc- The final state charmed hadron admixture is obtained from neutrino 
data from an emulsion experiment that has been corrected for improved knowledge of charmed hadron lifetimes 
p| . Extraction of CKM matrix elements from the data requires an independent determination of B^, this analysis 
uses the Particle Data Group charm semi-leptonic branching fractions convolved with the species production 
cross-sections just mentioned. The result is B^^' — 0.093 ± 0.009. 

In next-to-leading-order (NLO) QCD, additional diagrams complicate the expression for CTc and also spoil the 
factorization of the various components of cr^^. The charm production cross-section also becomes dependent on the 
QCD factorization and re-normalization schemes and their respective scales. 



III. EXPERIMENTAL MEASUREMENT AND ANALYSIS TECHNIQUE 



In an ideal situation one would like to present direct measurements of the differential charm production cross- 
sections doc'^'^^^ /dxdy at several different neutrino energies. Neither NuTeV nor its predecessor CCFR measures 
charm, but rather dimuons. The charm cross-section is thus related to the data by model-dependent corrections for 
charm fragmentation and decay and by experimental effects of resolution, acceptance, and neutrino flux. One way of 
handling these issues is to fit a parametric model directly to the data and extract parameters from the model. This 
approach was used in the past for LO QCD (l4|-|l| and NLO QCD |l| in the variable flavor ACOT |] scheme. 



The approach taken here begins with the same idea, a LO QCD parametric fit based on Eq. 0. Events passing 
selection criteria detailed in the next section are binned separately in and mode in the quantities 

Evis = Ef^i + Ef^2 + Eh AD , (3) 
^^^^ - 2Mii2 + EHAu) '"^'^ 

^^'s = p , p , (5) 

where E^i is the energy of the primary muon with the same lepton number as the beam, £"^2 is the energy of the 
other muon, E^ad is the observed hadronic energy in the calorimeter, and is the scattering angle of the primary 
muon. The "V^/S"' subscript indicates that these quantities differ from the true values of x, E, and z, due to the 
energy carried away by the neutrino from charm decay and due to detector smearing. Other quantities of interest for 
comparison purposes are 

Ef_i2 + Ehad , ,„s 

WIS = p , p , p and (6) 

Qvjs = 2MEvisxvisyvis- (7) 

A binned likelihood fit is performed which compares the data to a model composed of a charm source described by 
Eqs. 1^ and |^. The charm events are augmented with a contribution from dimuon production through tt/K decay in 
the charged-current neutrino interaction's hadron shower and then processed through a detailed Monte Carlo (MC) 
simulation of the detector and the same event reconstruction software used for the data. The MC dimuon sample 
is normalized to the data through use of the inclusive single muon event rates in v^^ and Df^ mode. The fit varies a 
common charm mass rric, branching fraction Be, and fragmentation parameter e for both modes, and two parameters 
for each mode, {k^, a,y) and {np, ap), that describe the magnitude and shape of the s and s quark PDFs. The strange 
sea parameters are defined by 

/ n2\ u{x,Q'^) +d{x,Q^) 
s[x,Q ) = K„ — — — [1 — X) and (8) 

-r n2\ u{x,Q^) +d{x,Q^) 

In these parameterization^ values of = Kp = 1 and = ap = would imply an SU(3)-flavor symmetric sea; 
previous measurements have yielded k values around 0.4, and a values consistent with zero (within large errors). The 



*This parameterization differs slightly from previous LO analyses in the definition of k,^ and Ka in the general case of 
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fragmentation process is described using the CoUins-Spiller fragmentation function [Q : 

e)^N [(1 - z)/z + ec{2 - z)/(l - z)] (1 + z^l - (l/z) - ec/{l - z)]-\ (10) 

The Peterson function |1^ was also tried but produced worse agreement between MC and data. 

The analysis proceeds based on the observation that the dimuon data are well described by the LO fit. The 
information on charm production from the LO fit is, in fact, not of great importance; our goal is to use the fit 
parameters to help construct a cross-section. This task is performed by forming a new grid in xvis, Vvis, a-nd Eyis 
in the data and MC, and a corresponding grid of x,y, and E in the MC. The dimuon cross-section is computed 
at the weighted center of each {x,y,E) bin i. The MC also predicts, with the result of the LO fit, the number of 
events in each {xviSy Uvis^ Eyjs) bin j. The MC can further be used to establish a correspondence between (x, y, E) 
and {xvis, Uvis, Eyis) bins; this is accomplished by finding the bin ji in {xvis, Uvis, Eyis) space that receives the 
largest fraction of events produced in {x, y, E) bin i. After this procedure the cross-section for (x, y, E) bin i is then 
determined by 



dxdy }{x.^v,E)^ KLOfit Je^2>e^2^,^ I dx dy dz )(x.,y,E)-wm (n) 

In this expression the left hand side represents the measured cross-section for dimuon production as a function of 
X, y, and E in bin i with the requirement that the second muon in the event exceeds the threshold used in the 
experiment, E^2 > E^2m\n- On the right hand side, Di is the number of data corresponding to bin i, and A/'/lq fit 
the corresponding number of events predicted by the LO fit. In the integrand, d^a^^{v^N — > fj,~ ^^X)/dx dy dz is 
taken from Eq. ^, and the integral over the fragmentation variable z and charmed hadron decay variable Q maintains 
the condition £'^2 > £'/j2min- The procedure for defining Di and A/^lq fit detailed further in Appendix ^ 

The end result is two tables, one each for z/^ and mode, of the "forward" dimuon cross-section (which is closest 
to what the experiment actually measures) . It will be shown later that these tables can be used to re-extract the LO 
fit parameters and can be combined with similar tables from the CCFR experiment. More details on the apparatus, 
event selection, and analysis procedure will be given first. 



IV. NUTEV EXPERIMENTAL DETAILS 



A. Detector and beamline 



The NuTeV (Fermilab-E815) neutrino experiment collected data during 1996-97 with the refurbished Lab E neutrino 
detector and a newly installed Sign-Selected Quadrupole Train(SSQT) neutrino beamline. Figure || illustrates the 
sign-selection optics employed by the SSQT to pick the charge of secondary pions and kaons which determine whether 
i^fi or i/p are predominantly produced. The SSQT produced beam impurities of (f^) events in i/^ mode (i/^ mode) 
at the 10-3 level. During NuTeV's run the primary production target received 1.13 x 10^^ and 1.41 x 10^* protons 
on target in neutrino and anti-neutrino modes, respectively, resulting in inclusive CC samples of 1.3 x 10® events in 
neutrino mode and 0.46 x 10® 9 inclusive CC events. The very low "wrong-flavor" backgrounds [ pO| imply that only 
one muon charge measurement is needed to make the correct assignments for Ef^i and £'^2 described above. 

The Lab E detector, described in detail elsewhere [2l| ], consists of two major parts, a target calorimeter and an iron 
toroid spectrometer. The target calorimeter contains 690 tons of steel sampled at 10 cm(Fe) intervals by 84 3m x 
3m scintillator counters and at 20 cm(Fe) intervals by 42 3m x 3m drift chambers. The toroid spectrometer consists 
of four stations of drift chambers separated by three iron toroid magnets that provide a px kick of 2.4 GeV/c. The 
toroid magnets were set to always focus the muon with same lepton number as the beam neutrino. Precision hadron 
and muon calibration beams monitored the calorimeter and spectrometer performance throughout the course of data 
taking. The calorimeter achieves a sampling-dominated resolution of oe/E — 2.4% © 87%/ ^/E and an absolute 
scale uncertainty of SE/E = 0.4%. The spectrometer's multiple-Coulomb-scattering-dominated muon momentum 
resolution is ap/p = 11%, and the muon momentum scale is known to Sp/p = 1.0%. 



and Oi? 7^ 0. The motivation for not using the older definitions (of the form s (x, Q^) ~ ^ {''^ (^1 Q^) + ^^)) " ^)"'' 
X „^ " z ) is to avoid a procedure that requires information about the PDF outside the experimentally 

dx(u(x,Q2)+d{x,Q^)){l-x)'''' 

accessible x range of the experiment. 
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B. Data selection 



A typical dimuon event has the characteristics shown in Fig. || In this figure, the toroid can be seen to focus the 
leading muon originating from the leptonic vertex and to de-focus the secondary muon, which originates most probably 
from charm decay. In the event shown both muons pass through the toroid, and both their signs are measured. In 
events where the sign of one muon is not measured, it is assumed to be opposite the one measured. Since the sign of 
the primary muon is known because of the sign selection of the SSQT, the measurement of the sign of only one muon 
is sufficient to identify the primary and secondary muon in the event. The rate of the same sign dimuon events with 
both muons toroid analyzed is very low pot] . 

Candidate opposite sign dimuon events, in both data and Monte Carlo, were selected using the following criteria: 



The event must occur in coincidence with the beam and fire the penetration trigger (charged-current interaction 
trigger) . 

The incident neutrino energy, E^, must be greater than 20 GcV and the energy of the hadronic shower, Ehad, 
greater than 10 GeV. 

In order to ensure event containment, only events occurring within an active fiducial volume are accepted: the 



transverse vertex positions {Vx,Vy) must be satisfy —127 cm< Vx^y < 127 cm and -^JV^ +Vy < 152.4 cm, and 

the longitudinal vertex position must lie between counters 15 and 80, which corresponds to 2.7 and 13.3 hadronic 
interaction lengths from the upstream and downstream ends of the calorimeter, respectively. 

• Two muons must be identified in the event and satisfy the following criteria: 

— The energy of each of the two muons, E^i,f^2, must be greater than 5 GeV. 

— The time obtained from fitting each track must be within 36 ns of the trigger time. 

— One of the muons must be toroid-analyzed and the energy of the toroid-analyzed muon at the entrance of 
the toroid must be greater than 3 GeV. 

— At least one toroid-analyzed muon must pass through at least 2/3 of the toroid. 

— The toroid-analyzed muons must hit the front face of the toroid inside a circle of radius Rff < 152 cm, 
and more than 80% of the path length of the muon must be in the toroid steel. 

• Finally, in order to remove mis-reconstructed events, a requirement is imposed on the reconstructed xyis 
kinematic variable: < xyis ^ 1- 

The final event sample contains 5102 z^^-induced and 1458 i^^-induced events. Of these, 2280/655 in v^Jv^^ mode 
have both muons reconstructed in the toroid spectrometer. All other events have only one. The mean Eyis of the 
events is 157.8 GeV, the mean Qyjg — 21.1 GeV^, and the mean xyis = 0.14. The overall reconstruction efficiency, 
including the detector acceptance, is ^ 60% for events with E^2 ~ 5 GeV, and ^ 80% when E^2 is above 30 GeV. 



C. Detector Simulation 

A hit-level Monte Carlo simulation of the detector based on the GEANT package was used to model the 
detector response and provide an accurate representation of the detector geometry. The Monte Carlo events were 
analyzed using the same reconstruction software used in the data analysis. The detector response in the simulation 
was tuned to both hadron and muon test beam data at various energies. To ensure the accurate modeling of the muon 
reconstruction efficiency the drift chamber efficiencies implemented in the simulation were measured in the data as a 
function of time. 

The prim ary neutrino interactions were generated using the LO QCD model and fragmentation function described 



in Sec. Electroweak radiative corrections based on the model by Bardin |23| were applied to this cross-section. 
The main background source from ordinary CC interactions in which a pion or kaon produced in the hadronic shower 
decays muonically was generated following a parameterization of hadron test beam muoproduction data for simulating 
secondary decays, and the LEPTO ||2^ package for the decays of primary hadrons 
produce such muons with momentum greater than 4 GeV/c is ~ 2 x 10~^ for events witl 



|25|]. The total probabihty to 

: 1 Ehad ~ 30 GeV, and « lO'^ 
for Ehad ^ 100 GeV. The contribution from the primary hadron decays is roughly two times larger than that from 
the secondary decays. 
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D. Neutrino Flux and Normalization 



The total flux, energy spectra, and composition for both neutrino and anti-neutrino beams are calculated using 
a Monte Carlo simulation of the beamline based on the DECAY TURTLE program ||2^ and production data from 
Atherton as parameterized by Malensek for thick targets. This flux is used to generate an inclusive charged- 
current interaction Monte Carlo sample using the GEANT based hit-level detector Monte Carlo described in section 
IV. C. 

The predicted flux is then tuned so that the inclusive charged-current interaction spectra in the Monte Carlo match 
the data. Selection criteria for this sample of inclusive charged-current interactions are exactly the same as those used 
to select the dimuon sample with the requirement for two muons removed. Flux corrections of up to 15% are applied 
in bins of neutrino energy and transverse vertex position to force the single muon data and Monte Carlo to agree. 
In addition, an overall factor is determined for each beam (neutrino or anti-neutrino) that absolutely normalizes the 
single muon Monte Carlo to the data. The dimuon Monte Carlo uses the flux determined with the above procedure; 
and it is absolutely normalized to the inclusive single muon charged-current data through the flux tuning procedure. 

The procedure used to tune the flux to the observed single muon rate is iterative since the event rate observed in the 
detector depends on the convolution of the cross-section with neutrino flux, and the result of the charm measurement 
has a small effect on the total cross-section. Corrections found from the single muon Monte Carlo/data comparisons 
are used in the dimuon Monte Carlo that is used to determine the dimuon cross-section, and thus the charm production 
cross-section within our LO model. The charm cross-section results are then used in the single muon Monte Carlo and 
the procedure is repeated until the flux parameters do not change (in practice the convergence is very fast). Figure ^ 
shows a comparison between data and Monte Carlo for the single muon (flux) sample; the level of agreement is very 
good. Note that muon and hadron energies are not adjusted separately in the flux tuning procedure. 

V. RESULTS 

A. NuTeV and CCFR Leading Order QCD fits 

The LO QCD fits were performed using three different parton distribution hmction (PDF) sets with their cor- 
responding QCD evolution kernels: GRV94LO Q and CTEQ4L0 |Q , as implemented in the PDF compilation 
PDFLIB [^, and a Buras-Gaemers parameterization ||3^] (BGPAR) that has been used extensively in this exper- 
iment and its CCFR predecessor. In the BGPAR case, an explicit Callan-Gross relation violation is implemented 

by replacing the term 1 - ^^^-^ in Eq. | with (1 + Rl){1 + (^f)~^{l -y- ^) + ^, where Rl, the ratio of 
longitudinal to transverse W^N cross-sections, is taken from a fit to electroproduction data Q. Table | hsts fit 
results with the rightmost column showing the combined for v and v modes. All three models have the same good 
level of agreement with the dimuon data. Figure |6| illustrates the quality of the BGPAR fit by comparing it to the 
data for the kinematic variables used directly in the fit; figure |^ shows a comparison for variables not used directly 
in the fit. One observes that differences in the choice of PDF parameterization result in different charm production 
parameters, indicating significant model dependence at LO. Because of this dependence, the most relevant quantity 
to extract is the dimuon production cross-section. 

B. CCFR Leading Order QCD fits 

The CCFR dimuon data set used in the CCFR LO and NLO analysis is used to extract LO QCD parameters 
using the same procedure described in the previous section. Results, shown in Table |l|, are consistent with the NuTeV 
fits except for the fragmentation function shape parameter e. This difference is caused by the fact that in NuTeV a 
much higher percentage of low energy muons is used, thus the zvis distribution has significant shape differences to 
that of CCFR. In addition, the uncertainty in the tt/K background parameterization is larger for low muon energies; 
this effect is reflected in the size of the error on the determination of e. 

Table || shows results of a combined fit to the data from both NuTeV and CCFR performed using the same 
procedure and the BGPAR PDF set. Since the fragmentation shape parameter is different for the two experiments, 
the Monte Carlo sample for each is reweighted to the appropriate e and then kept fixed in the combined fit to simplify 
the fitting procedure. 
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C. Comparing CCFR with NuTeV 



Different BGPAR PDF sets are used to analyze NuTeV and CCFR data. Thus, since we have shown that different 
PDF sets can produce different LO parameters, the numbers presented in Table || should not be compared directly 
with each other. Furthermore, the combined fit uses the NuTeV BGPAR PDF set, an approach that is not completely 
rigorous. Parameters in Table || should be compared to their corresponding values in Table [II, which gives the result 
of the fit to the dimuon cross-section table discussed in the next section. In this case, the use of any PDF set or model 
is equally appropriate, since the purpose is to compare the results of the fit to the dimuon data to those of the fit to 
the extracted cross-section tables. 

We caution readers against taking k and a parameters for strange seas extracted using the BGPAR sets and using 
them to construct s {x, Q^) /s {x, Q^) PDFs according to Eqs. ||-^ from a different PDF set. Since the NuTeV/CCFR 
PDFs are not publicly available, a safer course would be to take the strange sea parameters for CTEQ or GRV. Even 
then, it is important to match k and a from our fits with the appropriate PDF set. 



D. Systematic uncertainties 



Although the main thrust of the analysis is the extraction of the dimuon production cross-section, the various 
sources of systematic uncertainty are presented by listing their contributions to the LO fit parameters. This is done 
for reasons of clarity, since the individual systematic uncertainty contributions add too many entries in the cross- 
section tables. These methods are completely equivalent since the systematic uncertainty on the parameters of the 
LO fits propagates directly to the cross-section measurement. 

The main sources of the systematic uncertainties arise from modeling uncertainties in the Monte Carlo. The most 
significant are the tt/K decay background simulation, the detector calibration from the analysis of test beam hadron 
and muon data as a functions of energy and position, and the overall normalization. In addition, in the case of the 
BGPAR fits, the uncertainty on the longitudinal structure function is important. 

The systematic error sources are given in Table IV. For the combined NuTeV-|-CCFR fit, systematic errors begin 
to dominate the uncertainties for several fit parameters. 



E. NuTeV Cross-Section Tables 



As was shown in the previous section, charm cross-section parameters depend on the details of the charm production 
model used in the fit. In addition, a fragmentation and decay model must be used to extract the charm production 
cross-section from the observed dimuon rate, introducing further model dependence. These model dependences are 
exacerbated by the experimental smearing correction, due to the missing neutrino energy of the charmed hadron decay 
and the detector resolution, and the substantial acceptance correction for low energy muons. In order to minimize 
model dependencies, we choose to present our result in the form of a dimuon production cross-section. 

We have shown in the previous section that we can obtain very good description of our dimuon data, independent 
of the charm production model assumptions in our Monte Carlo. We thus limit the use of the Monte Carlo to correct 
experimental effects to the measured dimuon rate, and we limit the measurement of this rate to regions of phase 
space where the acceptance of the experiment is high. In this case, the only model dependence comes from potential 
smearing effects close to our acceptance cuts, i.e., the uncertainty on production of events produced with kinematic 
variables outside our cuts that smeared to reconstructed values within the cuts. The prediction for this kind of 
smearing depends on the underlying physics model, which is not well constrained by our data, since it involves phase 
space not accessed by our data. However, this model dependence is a second order effect. 

The dimuon cross-section extraction procedure depends on the ability of our Monte Carlo to describe the data. 
The "true" three-dimensional phase space (x, y, E) is divided into a number of bins. The grid is set up so that there 
is the same number of Monte Carlo events in each bin for any projection into one dimension. The cross-section for 
dimuon production with £^^2 ^ 5 GeV is calculated according to Eq. m at the center of gravity of each bin. The 
size of the bins in "visible" phase space (xy/s, yvis, Eyis) is defined in such a way so that there is a correspondence 
between a visible bin ji and a generated i established by the condition that bin ji contains at least 60% of the events 
from bin i. The relative error on the cross-section is taken to be equal to dg-, /(J2u = (^2^,/ \/ D (ji). Cross-section 
tables were produced using all three BGPAR, GRV, and CTEQ models; figure demonstrate the insensitivity of 
the cross-section to PDF choice. Further details of the binning procedure are defined in Appendix 
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F. Charm Production Fits Using the Dimuon Cross-Section Tables 



Using the cross-section tables involves similar steps as used in the direct data analysis, but with all detector and 
flux dependent effects removed. One must provide a model for charm production, fragmentation, and decay; construct 
the dimuon cross-section number for each entry in the table; and perform a ^ fit. The function should employ 
the statistical and systematic errors in each bin added in quadrature. 

One must also account for correlations between the various table entries. These correlations derive from our use 
of the LO fit to parameterize the data and from our method of binning; they are an inherent consequence of the 
incomplete kinematic reconstruction of the dimuon final state. We have adopted a pragmatic approach towards 
handling this issue. Rather than compute a large correlation matrix, we inflate the cross-section errors in each bin 
by a factor that is typically 1.4. This factor is chosen so that an uncorrelated fit to the tables returns the same 
parameter errors as a direct flt to the data. A consequence of this error inflation is an apparent best flt per bin 
(x^/bin) that is approximately 1/4 rather than 1. This low value results from overcounting the number of degrees of 
freedom (DOF) present in the table and is compensated for by using an effective DOF per bin determined via MC 
calculation and given along with cross section data in the tables. Further details of the determination of the effective 
DOF are presented in Appendix 

We tested this fitting procedure on the tables using the same BGPAR, fragmentation, and decay models used to 



obtain the table; Table III summarizes this study. Both the parameter values and their uncertainties obtained from 
fitting to the cross-section table agree with the corresponding values obtained by fitting directly to the data. It has 
also been verified that GRV94 parameters, for example, can be extracted from a cross-section table constructed with 
either the CTEQ or BGPAR model so as to agree with parameters obtained by fitting directly to the data. 

While our cross-checks in fitting the cross-section tables entail using the same physics model used in generating 
the tables, we emphasize that the table presents a set of physical observables which may be used to test any dimuon 
production model. For the most interesting case of dimuon production through charm, a typical model test would 
consist of the following steps^: 

• For each bin, generate a (large) ensemble of Nqen events with x, y, and Pc, where pc denotes the lab mo- 
mentum of the produced charm quark, according to the model charm production differential cross-section 
dcTc-mod/dx dy dp^ 

• Fragment the charm quarks into hadrons and decay the charmed hadrons (using, for example, PYTHIA), and 
determine NpAss the number of events which have a charm decay muon with E^2 > 5 GeV. 

• The cross-section to compare to the table value is then da^^^^^o^/dx dy = {NpAss /Nqen) x 
/ dpc doc-mod/dx dy dpc- 



• A fit should then be performed to minimize 



(da^^^n,od/dxdy-da+^^^^^/dxdy) 
^ = ^ ^ ^^-^^ (12) 

table-bins " stat ^ ^ syst 

in each beam mode with respect to the desired parameters in dac~mod/dx dy dpc- 

• The confidence level for the fit may be obtained by comparing the to the sum of effective DOF for table bins 
used in the fit. 

To provide further experimental information for model testing, the following kinematic quantities are given along 
with the cross-section in each bin: {Eh ad), the mean visible hadronic energy; (£'^2), the mean energy of the 
secondary muon; (p^2m)i mean square of the secondary muon's transverse momentum in the event scattering 
plane; and(p^2oiit)' the mean square of the secondary muons transverse momentum perpendicular to the event 
scattering plane. These quantit ies ar e computed from the dimuon data, with the LO fit used only for acceptance and 
smearing corrections. Tables V XVI contain the measurements^. Appendix ^contains a supplementary discussion of 
the cross sections at high x. 



tA simple PYTHIA implementation of the first three steps may be obtained at www-e815.fnal.gov, or by contacting the 
authors. 

■''These data may also be obtained electronically at www-e815.fnal.gov. 
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VI. SUMMARY 



We present a measurement of the dimuon production cross-section from an analysis of the data of the NuTeV 
neutrino DIS experiment at the Tevatron. NuTeV data are combined with an earher measurement from the CCFR 
experiment that used the same detector but a different beamhne. A leading order QCD analysis of charm production 
performed on the combined data yields the smallest errors to date on model parameters describing the charm mass, 
the size and shape of the strange sea, and the mean semi-muonic branching fraction of charm. The leading order QCD 
model describes NuTeV and CCFR data very well, but cross-section model parameters extracted vary depending on 
the particular choice of model. The extracted dimuon production cross-section, by contrast, is insensitive to the choice 
of leading order QCD model and should be of most use to the community of phenomenologists. 



APPENDIX A: CROSS-SECTION TABLE BINNING PROCEDURE 



In this analysis we report the differential cross-sections da [v^^/v^^Fe — s- /x^/i^X) / dxdy for forward secondary muons 
tabulated in bins of neutrino energy Bjorken scaling variable x, and inelasticity y. The measurement is obtained 
by using a LO Monte Carlo fit to the data to find the correspondence (mapping) between the "true" (unsmeared) 



and reconstructed (smeared) phase space, as discussed in Sec. This procedure maps the statistical fluctuations 
of the observed event rate to the "true" phase space bins where the cross-section is reported. The consistency of 
the procedure is checked by comparing fits of various models to the extracted cross-section tables to fits of the same 
models performed directly to the data. The criteria for the check is that the obtained central values and the errors 
on the model parameters are the same in both cases. In order to meet these criteria the binning of the smeared and 
the unsmeared phase space has to be appropriately selected. 

Usually in cross-section measurements, the procedure followed is to bin both the "true" and the reconstructed 
variables using the same grid, and select the bin size empirically so that for each bin the purity is maximized and 
the smearing contribution from other bins is minimized^ In such a method the correspondence between smeared and 
unsmeared bins is trivial since the same grid of bins is used. In the case of the dimuon cross-section measurement, a 
complication arises from the large smearing effects due to the missing neutrino energy in the reconstructed dimuon 
final state. Unlike detector resolution effects, this smearing is not a symmetric function of the "true," unsmeared 
variables, but rather an asymmetric mapping similar to electroweak radiative corrections. Because of this effect, we 
have followed a more elaborate procedure to map the visible phase space bins to those of the "true" phase space. This 
procedure allows us to obtain the desired high purity for each bin and also achieve stability of the result independently 
of the binning choice. In addition to mapping, our procedure allows us to take into account the significant bin-to-bin 
correlations which arise from the large smearing corrections without having to construct a full error matrix. This 
correlation matrix can be calculated, but it is too unwieldy to be useful; and it is difficult to incorporate effects of 
correlated systematic errors in a meaningful way. In our treatment, we estimate an effective number of degrees of 
freedom which allows us to obtain from an uncorrelated fit to the cross-section tables the same fit parameter errors 
as the ones obtained by directly fitting to the data, using the usual Ax'^ = 1 definition. 

We begin the description of the technique by defining more precisely the factors Di and A/^ lq fit from Eq. |ll|. In 
this equation, the expression on the left hand side represents the measured cross-section for dimuon production as a 
function of x, y, and E in bin i of the true phase space (with the requirement E^2 > Ef^2min)- On the right hand side, 
Di is the number of data events in bin i (which has to be determined using the Monte Carlo mapping procedure) , and 
A/7lo fit is the number of events predicted in bin i by the LO fit. To estimate the number of data events associated 
with bin i, we start by selecting bins in generated phase space {x, y, E)i, requiring equal number of Monte Carlo events 
in each projection of the true (generated) phase space, so the number of events in each (x, y, E)i bin is approximately 
the same. The visible phase space is divided using the same algorithm but with a much finer grid than the generated 
space. The mapping matrix M.ij makes the correspondence between visible and generated phase spaces: 

where Afij is the number of Monte Carlo events generated in bin i which end up in visible bin j and A/} is the total 



^Purity is defined here as the fraction of events which have both their unsmeared and smeared variables belonging to the same 
bin. 
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number of Monte Carlo events in visible bin j. The coverage fraction C in visible space is defined as 



N{C) 

C=Y. -^u/M, (Al) 

3 

where Mi is the number of Monte Carlo events in generated bin i. In the case of C = 1 the sum is performed over all 
visible bins; otherwise, the summation goes over the bins with the highest Mij/Mi ratios until the desired fractional 
coverage is obtained. Using the above definitions, and for a given coverage C, we can define the number of data events 
which belong to a given true phase space bin i (where the cross-section is reported) as 

D,=J2M^rDj. (A2) 

The number of Monte Carlo events in this generated bin, for a given coverage C, is redefined as 

A/^' = ^M,AG-- (A3) 

The cross-section error for each bin i should be proportional to the visible events "mapped" in that bin, so the 
following expression is used to assign it: 

6i = ^^m:Liju^/^^_ (A4) 

ax ay 

Note that the Monte Carlo statistics contribution comes from the total number of Monte Carlo events generated in 
bin i. The multiplicative factor Di/M^ in Eq. |ri| cancels out to first order the model dependence of the extracted 
cross-section and approximately transfers the statistical fluctuation in the visible phase space to the true phase space. 

The procedure as described above is incomplete, since the true bins i are not statistically independent. As we 
stated in the introduction to this section, we do not calculate the full error matrix, but rather estimate an effective 
(independent) number of degrees of freedom. It is possible to estimate the number of independent degrees of freedom 
by calculating the contribution to the total number of degrees of freedom from each bin 

DOF. = (A5) 
It is obvious that the effective number of degrees of freedom depends on the selected coverage fraction C, and should 



decrease as C increases. Figure 12 shows the number of effective degrees of freedom (DOF) as a function of the 
coverage area C (solid curve). The dotted curve in Fig. shows the obtained as a result of the fit to the table. We 
conclude that our method produces the correct number for / DOF, if the effective number of degrees of freedom is 
used, for a coverage fraction between 55% and 90%. 

The coverage area percentage used for our reported cross-section result is based on a Monte Carlo study. In this 
study a Monte Carlo sample is used to produce "cross-section" tables and then fits are performed to the tables and 
directly to the sample. Using as guidelines the criteria that there should be no pull on fit parameters in a fit to the 
cross-section table versus a direct fit, and that an uncorrelated fit to the tables should yield the same fit parameter 
errors as the ones obtained by a direct fit, we selected a coverage area C = 60%. The effective number of degrees of 
freedom which corresponds to this value should be used with all fits performed on the cross-section tables we present 
in this article. 



APPENDIX B: THE HIGH X REGION 



This section presents a supplementary investigation of the high x region [x > 0.5). The objective of this study is to 
ascertain whether there is any indication of an enhancement in the cross-section that we may be missing due to our 
use of wide x bins (dictated by the low observed event rate at high x). Such an enhancement could be caused by an 
unusually large strange sea, particularly in neutrino mode, which has been advocated to resolve certain discrepancies 
between inclusive charged lepton and neutrino scattering | |3^ . Previous dimuon analyses p^ , p^ may have missed this 
effect due to dependence on the particular model used to parameterize the strange sea distribution. 

In order to minimize model-dependent corrections, we report our high x cross-section measurements as fractions of 
the total dimuon cross-section. Similarly, to quote a limit for the a; > 0.5 cross-section we use the observed data rate 
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for Xvis > 0.5. This is a conservative way to set a limit, since by the kinematic effect of the missing decay neutrino 
energy, the contribution to a given xyis bin always comes from x < xyis- 

The cross-section ratio of the dimuon cross-section for x > 0.5 to the total dimuon cross-section in a given energy 
bin can be expressed as 



2m 
^a;>0.5 
2ii 



N. 



x>0.5 



M 



a 



Ntot Me>5' 



(Bl) 



where Nx>o.5 is the number of observed events for xvis > 0.5, Ntot is the total number of observed dimuon events, 
M is the Monte Carlo prediction with all experimental cuts applied, and Me>5 is the Monte Carlo prediction without 
the tt/K decay contribution and wit h onl y the -E^2 > 5 GeV cut applied. For this study, we use the same data 
selection criteria described in section IVB, exept for the xyis selection, which is changed to: < xyis ^ 2. 



In the anti-neutrino dimuon data sample we define three energy bins and we record the number of the observed 
dimuon events with xvis > 0.5 in the data, together with the Monte Carlo prediction and all the relevant information 
for Eq. Bl, Here the Monte Carlo prediction is very well constrained by our dimuon data in the full x range, since the 
observed rate is mostly due to scattering on s quarks. The results are presented in Table pCVII . The systematic and 
statistical error on the Monte Carlo prediction is negligible for this discussion. Treating the Monte Carlo prediction 



as a "background," and using Eq. p31| , we set cross-section ratio upper limits at 90% CL, for any additional s ource 
of a; > 0.5 dimuons, of 0.0012, 0.007, and 0.009, respectively, in each one of the energy bins defined in Table XVII 
(counting from lower to higher energy bin). 

We follow the same procedure in the neutrino data sample. Here, there is an additional complication in the 
interpretation of the result since a large contribution from valence quark events is expected. The Monte Carlo 
prediction for the rates of the non-strange sea component is not directly constraine d by ou r dimuon data, but rather 
by inclusive structure function measurements. The results are presented in Table XVIII. It is worth noticing that 
within the model and PDF sets used in this analysis (BGPAR) the Monte Carlo prediction for the contribution of 
the strange sea is only on the order of 3.7% of the total rate for xvis > 0.5; most of the contribution (69%) in this 
model comes fro m the valence quarks. 

Using Eq. Bl and treating the Monte Carlo prediction as a "background" source, we set 90% CL li mits fo r an 
additional cross-section source at x > 0.5. We find that for the first and last energy bins in table XVIII this 



additional source cannot be larger than 0.006 and 0.013 of the total dimuon cross-section, while for the 153.9-214.1 
bin there is less than 5% probability that there is an additional source consistent with our data (note that we have a 
1.75(7 negative yield compared to our background prediction). 
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model 


ruc (GcV/c") 


e 


Be (%) 


X'/DOF 


BGPAR 


1.33 ±0.19 ±0.10 


2.07 ±0.31 ±0.64 


11.40 ±1.08 ±1.15 


105/112 


GRV 


1.65 ±0.18 ±0.09 


2.09 ±0.31 ±0.64 


11.11 ± 1.51 ± 1.60 


101/112 


CTEQ 


1.63 ±0.17 ±0.09 


2.07 ±0.31 ±0.63 


10.70 ±1.66 ±1.76 


100/112 



model 


AC 


AC 


Q 


a 


BGPAR 


0.32 ± 0.06 ± 0.04 


0.37 ±0.05 ±0.04 


-1.10 ± 1.05 ±0.59 


-2.78 ±0.42 ±0.40 


GRV 


0.37 ±0.05 ±0.03 


0.37 ±0.06 ±0.06 


0.87 ±1.25 ±0.71 


0.28 ±0.44 ±0.42 


CTEQ 


0.44 ± 0.06 ± 0.04 


0.45 ±0.08 ±0.07 


1.17 ±1.20 ±0.68 


1.08 ±0.44 ±0.41 



TABLE I. Results of LO fits to NuTeV data. The first error is statistical, the second systematic. 



SET 


m, (CeV/c^) 


e 


(%) 


AC 


AC 


a 


a 


NuTeV 


1.33 ±0.19 


2.07 ±0.31 


11.40 ± 1.08 


0.32 ±0.06 


0.37 ±0.05 


-1.10 ± 1.05 


-2.78 ±0.42 


CCFR 


1.20 ±0.23 


0.88 ±0.12 


11.43 ±0.95 


0.31 ±0.05 


0.36 ±0.05 


3.14 ±0.91 


3.46 ±0.73 


combined 


1.38 ±0.13 




11.57 ±0.70 


0.35 ± 0.04 


0.41 ±0.04 


-0.77 ±0.66 


-2.04 ±0.36 



TABLE IL Results of LO fits to NuTeV, CCFR, and combined data set, using the BGPAR PDF set. Only the statistical 
errors are shown. 



SET 


m, (GeV/c^) 


B, (%) 


AC 


AC 


a 


a 


NuTcV 


1.30 ± 0.22 


10.22 ± 1.11 


0.38 ±0.08 


0.39 ±0.06 


-2.07 ±0.96 


-2.42 ±0.45 


CCFR 


1.56 ± 0.24 


12.08 ±0.99 


0.28 ±0.05 


0.33 ±0.05 


3.85 ± 1.17 


3.30 ±0.83 


NuTcV± 


1.40 ± 0.16 


11.00 ±0.71 


0.36 ±0.05 


0.38 ±0.04 


-1.12 ± 0.73 


-2.07 ±0.39 



TABLE III. Results of LO fits to the cross-section tables extracted from the NuTeV, CCFR, and combined data sets. 
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nic (^Lie V /c j 


e 


tic [/O) 


K 


K 


Oi 


(X 


v 7r//v (io/cj 


U.Uzz 


U.Oi 


U.Oi 


A ni Q 
U.Uio 


A AQ 1 

U.Uoi 


A A1 
U.UI 


A AC 


V TV / K (Zi70 ) 


U.UUD 


U.io 


U.UO 


A AAI 
U.UUi 


A A1 ^7 


A AT 
U.Ui 


A 1 ^ 
U.i / 




0.037 


0.09 


0.17 


0.001 


0.010 


0.48 


0.26 


fj, energy scale (1%) 


0.080 


0.33 


0.74 


0.036 


0.023 


0.25 


0.24 


Hadron energy scale (0.4%) 


0.012 


0.08 


0.02 


0.005 


0.003 


0.01 


0.04 


MC statistics 


0.047 


0.02 


0.31 


0.012 


0.006 


0.23 


0.01 


Flux 


0.010 


0.01 


0.07 


0.001 


0.000 


0.03 


0.03 


Systematic Error 


0.104 


0.64 


1.15 


0.043 


0.043 


0.59 


0.40 



TABLE IV. Systematic error sources for the LO-QCD fit to the NuTeV data. 



•7 

X 


X 


y 


d.c,l„ \- 


{I^iiad) 


(^V ) 




\l>\, ■„,.,) 


0.64 


0.021 


0.334 


0.419 ±0.071 ±0.003 


18.3 ± 1.0 


7.9 ±0.5 


0.40 ± 0.55 


0.14 ±0.24 


0.45 


0.058 


0.334 


0.538 ± 0.090 ± 0.022 


18.6 ± 1.2 


9.0 ±0.6 


1.08 ± 1.14 


0.76 ± 1.49 


0.44 


0.102 


0.334 


0.427 ± 0.069 ± 0.007 


17.9 ± 1.2 


8.3 ± 0.7 


1.19 ±0.78 


0.26 ± 0.48 


0.46 


0.168 


0.334 


0.323 ± 0.049 ± 0.008 


18.7 ±1.2 


8.2 ±0.5 


0.74 ± 0.28 


0.14 ±0.25 


0.67 


0.324 


0.334 


0.132 ± 0.019 ± 0.003 


19.4 ±0.9 


8.5 ± 0.4 


0.95 ± 0.09 


0.07 ± 0.04 


0.64 


0.021 


0.573 


0.774 ±0.107 ±0.015 


36.3 ± 1.4 


10.9 ±0.8 


0.19 ±0.09 


0.09 ±0.03 


0.46 


0.058 


0.573 


0.808 ± 0.108 ± 0.027 


34.3 ± 1.5 


10.9 ±0.7 


0.29 ±0.14 


0.11 ±0.07 


0.47 


0.102 


0.573 


0.792 ±0.103 ±0.012 


36.5 ± 1.5 


10.0 ±0.6 


1.08 ±0.97 


0.44 ± 0.96 


0.50 


0.168 


0.573 


0.471 ±0.060 ±0.017 


34.8 ± 1.6 


10.6 ±0.7 


0.85 ± 0.72 


0.25 ±0.43 


0.62 


0.324 


0.573 


0.198 ±0.027 ±0.003 


34.6 ± 1.6 


11.1 ±0.7 


0.74 ±0.11 


0.07 ±0.03 


0.58 


0.021 


0.790 


0.795 ±0.126 ±0.096 


49.6 ± 3.0 


14.6 ± 1.4 


0.26 ± 0.35 


0.10 ±0.05 


0.43 


0.058 


0.790 


0.894 ±0.133 ±0.029 


52.8 ± 2.5 


13.5 ± 1.5 


0.34 ± 0.35 


0.24 ±0.44 


0.41 


0.102 


0.790 


0.826 ±0.123 ±0.027 


52.1 ± 2.5 


11.2 ± 1.0 


0.24 ±0.21 


0.09 ±0.08 


0.52 


0.168 


0.790 


0.706 ±0.108 ±0.005 


52.9 ± 2.5 


12.6 ± 1.2 


0.30 ±0.15 


0.10 ±0.08 


0.58 


0.324 


0.790 


0.210 ±0.038 ±0.004 


49.0 ±2.9 


12.8 ± 1.2 


0.52 ±0.18 


0.09 ± 0.04 



TABLE V. NuTeV forward differential cross-section for Vij,N — > fj,~ iX^X at E ^ 90.18 GeV. The forward cross-section 
requires > 5 GeV, and the cross-section values should be multiplied by x G%MEj'K. The first error given for the 
cross-sections is statistical and the second systematic. Units are in GeV or GeV , where appropriate, for the averages of the 
kinematic quantities. 
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X 


y 


dxdy 


(Ehad) 


(e„+) 




<p2, . . ) 

\^T^+out/ 


0.66 


0.021 


0.334 


1.013 ± 0.148 ± 0.022 


39.1 ± 1.8 


12.1 ± 1.0 


0.57 ± 1.32 


0.24 ± 0.49 


0.45 


0.058 


0.334 


0.837 ± 0.127 ± 0.018 


39.6 ± 2.1 


12.2 ± 1.1 


0.33 ± 0.09 


0.08 ± 0.04 


0.40 


0.102 


0.334 


0.737 ± 0.111 ± 0.009 


39.8 ± 2.3 


12.2 ± 1.1 


0.57 ± 0.17 


0.08 ± 0.03 


0.39 


0.168 


0.334 


0.484 ± 0.071 ± 0.010 


38.5 ± 2.2 


11.3 ± 0.9 


0.69 ± 0.15 


0.11 ± 0.10 


0.57 


0.324 


0.334 


0.212 ± 0.028 ± 0.005 


39.5 ± 1.8 


12.7 ± 0.8 


1.44 ± 0.20 


0.08 ± 0.03 


0.57 


0.021 


0.573 


1.304 ± 0.196 ± 0.015 


69.0 ± 2.3 


18.1 ± 1.8 


0.19 ± 0.07 


0.11 ± 0.06 


0.43 


0.058 


0.573 


1.161 ± 0.176 ± 0.021 


71.0 ± 2.5 


16.7 ± 1.5 


0.25 ± 0.08 


0.08 ± 0.04 


0.35 


0.102 


0.573 


1.140 ± 0.178 ± 0.019 


72.9 ± 3.2 


19.2 ± 2.1 


0.63 ± 0.24 


0.26 ± 0.78 


0.42 


0.168 


0.573 


0.685 ± 0.107 ± 0.007 


75.0 ± 3.2 


14.8 ± 1.6 


1.34 ± 1.26 


0.15 ± 0.27 


0.52 


0.324 


0.573 


0.242 ± 0.038 ± 0.004 


76.8 ± 3.3 


13.6 ± 1.6 


0.74 ± 0.24 


0.07 ± 0.04 


0.43 


0.021 


0.790 


1.267 ±0.179 ±0.025 


101.8 ±3.6 


22.0 ±2.8 


0.20 ± 0.09 


0.11 ±0.05 


0.34 


0.058 


0.790 


1.301 ±0.176 ±0.031 


101.5 ±3.2 


21.1 ± 2.4 


0.29 ±0.11 


0.09 ± 0.04 


0.32 


0.102 


0.790 


1.072 ±0.153 ±0.023 


98.6 ±3.6 


21.9 ±3.1 


0.37 ±0.17 


0.09 ±0.05 


0.37 


0.168 


0.790 


0.788 ±0.118 ±0.020 


101.8 ±4.3 


20.4 ± 3.3 


0.72 ± 0.65 


0.09 ± 0.06 


0.41 


0.324 


0.790 


0.251 ± 0.040 ± 0.005 


100.4 ± 4.2 


20.2 ±3.1 


0.69 ± 0.40 


0.08 ± 0.07 



TABLE VI. NuTeV forward differential cross-section for z/^iV — > jX^X at i? ~ 174.37 GeV. The forward cross-section 

requires > 5 GeV, and the cross-section values should be multiplied by y^jj x G%AIE/tt. The first error given for the 

cross-sections is statistical and the second systematic. Units are in GeV or GeV , where appropriate, for the averages of the 
kinematic quantities. 





X 


y 


dxdy H" 


{Ehad) 




{PTi,+ in) 


{PTfi+out) 


0.79 


0.021 


0.334 


1.154 ±0.167 ±0.015 


54.8 ± 2.9 


15.1 ± 1.5 


0.82 ± 1.29 


0.09 ± 0.05 


0.54 


0.058 


0.334 


1.306 ±0.196 ±0.034 


59.3 ±3.1 


14.2 ± 1.3 


0.31 ±0.10 


0.29 ±0.40 


0.46 


0.102 


0.334 


0.913 ±0.132 ±0.017 


61.2 ±4.5 


13.2 ± 1.6 


0.40 ±0.11 


0.11 ±0.14 


0.41 


0.168 


0.334 


0.712 ±0.103 ±0.013 


65.7 ±3.3 


14.1 ±1.4 


0.82 ±0.18 


0.09 ± 0.03 


0.63 


0.324 


0.334 


0.278 ± 0.036 ± 0.003 


69.9 ±3.5 


13.8 ± 1.2 


1.29 ±0.30 


0.11 ±0.10 


0.68 


0.021 


0.573 


1.487 ± 0.234 ± 0.010 


97.7 ±3.8 


20.1 ±3.0 


0.19 ±0.09 


0.14 ±0.16 


0.47 


0.058 


0.573 


1.419 ±0.235 ±0.016 


105.7 ±3.9 


19.7 ±2.6 


0.83 ± 0.98 


0.10 ±0.09 


0.46 


0.102 


0.573 


1.018 ±0.166 ±0.011 


105.5 ±6.1 


21.2 ± 3.9 


0.84 ±0.62 


0.13 ±0.52 


0.42 


0.168 


0.573 


0.700 ±0.116 ±0.005 


108.1 ±4.9 


17.5 ±3.4 


0.78 ±0.34 


0.10 ±0.07 


0.60 


0.324 


0.573 


0.294 ± 0.046 ± 0.003 


109.6 ±4.5 


18.6 ± 2.7 


0.86 ±0.27 


0.08 ±0.04 


0.56 


0.021 


0.790 


1.656 ±0.222 ±0.060 


148.9 ± 4.9 


22.6 ± 3.3 


0.19 ±0.11 


0.16 ±0.20 


0.43 


0.058 


0.790 


1.546 ± 0.216 ± 0.013 


144.3 ± 5.5 


27.6 ± 4.5 


0.39 ± 0.20 


0.13 ±0.17 


0.40 


0.102 


0.790 


1.211 ±0.174 ±0.046 


149.8 ± 7.2 


25.4 ±4.4 


1.25 ± 1.19 


0.26 ±0.51 


0.36 


0.168 


0.790 


0.968 ±0.141 ±0.011 


150.3 ±6.9 


22.2 ±3.5 


0.52 ±0.32 


0.48 ± 1.98 


0.37 


0.324 


0.790 


0.311 ±0.047 ±0.005 


154.5 ± 6.3 


19.2 ±3.6 


0.50 ± 0.26 


0.07 ± 0.04 



TABLE VII. NuTeV forward differential cross-section for z/^iV — » fj/X at ~ 244.72 GeV. The forward cross-section 

requires E^+ > 5 GeV, and the cross-section values should be multiplied by y±j x G%ME /-k . The first error given for the 
cross-sections is statistical and the second systematic. Units are in GeV or GeV , where appropriate, for the averages of the 
kinematic quantities. 
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X 


y 


dxdy 


(Ehad) 


\ M / 




(pI 

\^Tfi out I 


0.61 


0.016 


0.356 


0.403 ± 0.101 ± 0.004 


16.9 ± 1.3 


7.5 ± 0.8 


0.11 ± 0.07 


0.06 ± 0.06 


0.39 


0.044 


0.356 


0.393 ± 0.101 ± 0.019 


15.9 ± 1.9 


8.2 ± 1.0 


1.34 ± 1.80 


0.59 ± 1.57 


0.34 


0.075 


0.356 


0.417 ± 0.107 ± 0.007 


16.5 ± 1.6 


7.9 ± 1.1 


0.42 ± 0.33 


0.11 ± 0.13 


0.37 


0.117 


0.356 


0.275 ± 0.067 ± 0.025 


16.5 ± 1.5 


8.4 ± 0.7 


0.37 ± 0.19 


0.12 ± 0.13 


0.63 


0.211 


0.356 


0.080 ± 0.018 ± 0.003 


18.5 ± 1.9 


8.0 ± 0.9 


0.74 ± 0.23 


0.08 ± 0.06 


0.61 


0.016 


0.586 


0.463 ± 0.097 ± 0.030 


30.6 ± 2.2 


13.0 ± 1.9 


0.19 ± 0.09 


0.05 ± 0.04 


0.39 


0.044 


0.586 


0.533 ± 0.110 ± 0.007 


31.2 ± 2.5 


10.3 ± 1.4 


0.17 ±0.08 


0.11 ± 0.05 


0.42 


0.075 


0.586 


0.621 ± 0.125 ± 0.026 


33.5 ± 2.2 


10.7 ± 1.4 


0.32 ± 0.17 


0.05 ± 0.05 


0.44 


0.117 


0.586 


0.357 ± 0.071 ± 0.008 


31.5 ± 2.4 


11.1 ± 1.3 


0.35 ± 0.13 


0.09 ± 0.06 


0.63 


0.211 


0.586 


0.183 ±0.036 ±0.003 


32.2 ± 2.3 


10.7 ±1.2 


0.93 ± 0.39 


0.07 ± 0.04 


0.53 


0.016 


0.788 


0.623 ± 0.161 ± 0.009 


43.7 ±3.8 


12.1 ± 1.6 


0.14 ±0.13 


0.13 ±0.20 


0.39 


0.044 


0.788 


0.623 ±0.154 ±0.016 


44.8 ± 4.3 


10.1 ± 1.8 


0.10 ±0.08 


0.04 ± 0.03 


0.42 


0.075 


0.788 


0.520 ±0.124 ±0.004 


43.4 ±3.7 


13.4 ±3.3 


0.19 ±0.11 


0.10 ±0.09 


0.41 


0.117 


0.788 


0.399 ± 0.098 ± 0.016 


41.1 ± 5.2 


11.9 ±2.7 


0.30 ±0.21 


0.09 ±0.10 


0.64 


0.211 


0.788 


0.153 ±0.039 ±0.004 


46.8 ± 5.1 


10.0 ± 2.4 


0.30 ± 0.26 


0.07 ± 0.06 



TABLE VIII. NuTeV forward differential cross-section for v^^N — > \x~ X at i? ~ 78.98 GeV. The forward cross-section 

requires _E^- > 5 GeV, and the cross-section values should be multiphed by G\ME j-K. The first error given for the cross-sections 
is statistical and the second systematic. Units are in GeV or GeV^, where appropriate, for the averages of the kinematic 
quantities. 





X 


y 


dxdy H" 


{Ehad) 




{PT/i-in) 


{PTii-out) 


0.68 


0.016 


0.356 


0.793 ±0.171 ±0.017 


32.8 ±2.5 


12.7 ±2.0 


0.50 ± 0.62 


0.78 ± 1.84 


0.42 


0.044 


0.356 


0.744 ±0.172 ±0.024 


36.1 ±3.3 


12.2 ± 2.0 


0.23 ± 0.12 


0.08 ±0.06 


0.37 


0.075 


0.356 


0.752 ±0.174 ±0.013 


33.6 ± 2.9 


11.2 ± 1.6 


0.33 ±0.11 


0.21 ±0.36 


0.38 


0.117 


0.356 


0.427 ± 0.097 ± 0.009 


36.7 ±3.7 


10.3 ± 1.4 


0.73 ± 0.27 


0.07 ± 0.06 


0.56 


0.211 


0.356 


0.140 ±0.030 ±0.003 


38.5 ±4.1 


8.8 ±0.8 


0.61 ± 0.22 


0.06 ±0.05 


0.53 


0.016 


0.586 


1.338 ±0.317 ±0.025 


58.5 ± 3.0 


16.4 ± 2.2 


0.09 ± 0.07 


0.07 ± 0.04 


0.37 


0.044 


0.586 


0.930 ± 0.223 ± 0.048 


66.3 ± 4.4 


13.2 ± 2.5 


0.12 ±0.09 


0.09 ±0.12 


0.36 


0.075 


0.586 


0.744 ±0.173 ±0.017 


60.1 ±4.5 


16.7 ±2.5 


0.27 ±0.13 


0.08 ± 0.08 


0.34 


0.117 


0.586 


0.653 ±0.159 ±0.005 


64.0 ± 3.9 


14.7 ±2.9 


0.38 ±0.18 


0.08 ± 0.07 


0.50 


0.211 


0.586 


0.257 ± 0.063 ± 0.006 


59.8 ±4.6 


14.4 ±3.2 


0.69 ± 0.32 


0.05 ±0.04 


0.46 


0.016 


0.788 


0.915 ±0.198 ±0.013 


83.0 ± 5.1 


18.1 ±3.4 


0.12 ±0.09 


0.09 ±0.14 


0.37 


0.044 


0.788 


1.106 ±0.240 ±0.042 


84.3 ± 6.4 


22.0 ±4.1 


0.30 ±0.18 


0.09 ± 0.07 


0.33 


0.075 


0.788 


0.775 ±0.172 ±0.041 


87.3 ± 5.8 


18.4 ±4.6 


0.23 ±0.31 


0.07 ±0.05 


0.36 


0.117 


0.788 


0.547 ±0.121 ±0.015 


88.2 ± 6.5 


15.7 ±4.6 


0.17 ±0.11 


0.06 ± 0.06 


0.44 


0.211 


0.788 


0.297 ± 0.075 ± 0.004 


78.7 ±5.6 


17.3 ± 3.6 


0.31 ±0.13 


0.07 ± 0.04 



TABLE IX. NuTeV forward differential cross-section for P^N —^ jj/ n~ X at S ~ 146.06 GeV. The forward cross-section 

requires E^- > 5 GeV, and the cross-section values should be multiphed by G%M E/tv. The first error given for the cross-sections 
is statistical and the second systematic. Units are in GeV or GeV^, where appropriate, for the averages of the kinematic 
quantities. 
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^2 

X 


X 


y 


dxdy 1^ 


\J^HAD} 




/n2 \ 


/r? \ 
\i^Tii- out/ 


n on 


O.Ui6 


0.006 


i.04D ± U.22o ± O.UiO 


00.0 ± 0.0 


14.0 ± 2.4 


Alt; _l_ A 1 A 

U.iD ± U.iU 


A 1 A ^ A AC 

0.10 ± U.Uo 


n 1 
U.Oi 


n n /I /I 
U.U44 


U.oOD 


i.ioo ± U.Z04 zh U.Ui / 


iA ^ A A 

00. U ± 4.4 


1 c n ^ A 
io.y ± Z.4 


A 97 _l_ A OQ 
U.O ( ± U.ZO 


A AQ ^ A A /I 
U.Uo ± U.U4 


n /1 1 
U.4i 


U.U/0 


0.000 


A c c _i_ nine _i_ n m A 
U.OOO ± U.i90 ± U.UiU 


C7 _1_ G A 

01 .2 ± D.4 


1 ^; _i_ c 
iD.2 ± 0.0 


A CO _1_ A 1 A 
U.OZ ± U.i9 


A A A _1_ A A'y 
U.U9 ± U.Vl 


U.44 


nil'? 


0.000 


A /lot? _1_ A An /I _1_ A AO A 

U.4ZD ± U.U94 ± U.UoU 


C7 1 _1_ '7 A 


1 Q C _1_ 

io.6 ± z.o 


A CO _l_ A o'y 
U.OO ± U.OI 


A A'? _1_ A AC 

U.U/ ± U.Uo 


n fin 
U.uU 


noil 
U.zii 


U.OOD 


A QQ1 _l_ A (AytA _l_ (A AOA 

u.ooi ± u.u/u ± u.uzy 


Di.y ± D.y 


ii.O ± Z.O 


1 /I _l_ A AO 

i.4D ± u.yz 


A AO -X- A A/1 
U.Uo ± U.U4 


U.DU 


U.UlO 


U.OOO 


1 /I CA _1_ A QOA _1_ A A1 

1.4oy ± U.ooU ± U.Ulz 


AC A _1_ £! T 
90. U ± D. / 


170 _1_ 

i I.Z ± 0.0 


A 1 /I _1_ All 

U.14 ± u.ll 


A A/I _1_ A A /I 

U.U4 ± U.U4 


U.4o 


U.U44 


U.OOD 


1 111 _1_ A 001 _1_ A A1 1 
l.iii ± U.zoi ± U.Ull 


AC A _1_ C A 

95. U ± 0.4 


1 n C _1_ Q A 

19.0 ± 0.9 


A 07 _1_ A 01 

U.zr ± U.zl 


A _1_ A AC 

U.Ul ± U.Uo 


n /1 1 
U.4i 


U.U(0 


U.OOD 


u.yyo ± u.zoo ± u.uo( 


CA _L 7 /I 

oy.z ± /.4 


OC 1 r n 
ZO.O ± O.Z 


A QA _l_ A 00 
U.Oy ± U.ZZ 


A 1 C ^ A 1 

U.lO ± U.lD 


U.4z 


U.ii ( 


U.056 


A 'yo'y _1_ A OAO _1_ A AT 1 

U.ioi ±U.zUz±U.Uii 


A'y _l_ 1 A 

y /.o ± io.o 


/I A _1_ 7 A 

24.y ± (.y 


1 /t _1_ 1 AG 

l.z4 ± l.lo 


A AO _1_ A 1 A 

U.Uo ± U.IU 


0.56 


0.211 


0.586 


0.303 ± 0.077 ± 0.007 


101.4 ±6.3 


18.6 ±5.8 


1.17 ±0.58 


0.16 ±0.33 


0.59 


0.016 


0.788 


1.125 ±0.243 ±0.018 


135.4 ± 12.9 


26.5 ±6.1 


0.15 ±0.13 


0.09 ± 0.09 


0.46 


0.044 


0.788 


1.433 ±0.295 ±0.012 


132.8 ± 10.2 


25.2 ± 6.0 


0.17 ±0.15 


0.08 ±0.11 


0.41 


0.075 


0.788 


1.258 ± 0.268 ± 0.033 


129.7 ±9.7 


26.3 ± 7.3 


0.52 ±0.54 


0.08 ±0.09 


0.44 


0.117 


0.788 


0.693 ±0.154 ±0.022 


134.2 ± 14.4 


25.8 ± 7.4 


0.70 ± 0.62 


0.06 ±0.07 


0.58 


0.211 


0.788 


0.219 ± 0.050 ± 0.010 


132.4 ± 18.0 


20.5 ± 5.3 


0.68 ± 0.35 


0.15 ±0.17 



TABLE X. NuTeV forward differential cross-section for i/^/V^ — > jx^ ix~ X at -B ~ 222.14 GeV. The forward cross-section 

requires E^- > 5 GoV, and the cross-section values should be multiplied by G'pM E/tv. The first error given for the cross-sections 
is statistical and the second systematic. Units are in GeV or GeV^, where appropriate, for the averages of the kinematic 
quantities. 





X 


y 


dxdy 


0.54 


0.023 


0.320 


0.589 ±0.129 ±0.018 


0.47 


0.057 


0.320 


0.619 ±0.110 ±0.004 


0.44 


0.100 


0.320 


0.472 ± 0.082 ± 0.004 


0.46 


0.167 


0.320 


0.494 ± 0.081 ± 0.003 


0.69 


0.336 


0.320 


0.159 ±0.024 ±0.001 


0.59 


0.023 


0.570 


1.367 ±0.223 ±0.019 


0.43 


0.057 


0.570 


1.017 ±0.161 ±0.008 


0.50 


0.100 


0.570 


0.706 ±0.103 ±0.023 


0.52 


0.167 


0.570 


0.472 ±0.067 ±0.010 


0.61 


0.336 


0.570 


0.226 ±0.034 ±0.003 


0.42 


0.023 


0.795 


1.406 ± 0.223 ± 0.047 


0.39 


0.057 


0.795 


1.361 ±0.203 ±0.016 


0.46 


0.100 


0.795 


1.138 ±0.172 ±0.014 


0.48 


0.167 


0.795 


0.812 ±0.136 ±0.018 


0.49 


0.336 


0.795 


0.211 ±0.040 ±0.003 



TABLE XI. CCFR E744/E770 forward differential cross-section for z/^iV jT ij+ X at E 109.46 GeV. The forward 

cross-section requires E^+ > 5 GeV, and the cross-section values should be multiplied by G%M E/tv. The first error given for 
the cross-sections is statistical and the second systematic. Units are in GeV or GeV^, where appropriate, for the averages of 
the kinematic quantities. 
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v2 

A. 


X 


V 


dxdy 1^ 


0.58 


0.023 


0.320 


1 319 + 231 + 008 


0.45 


0.057 


0.320 


1.161 ± 0.192 ± 0.036 


0.46 


0.100 


0.320 


qrin + n 1 47 + f) nnti 


0.41 


0.167 


0.320 


f) Kcy _(_ Q nnn _(_ n ririo 


0.68 


0.336 


0.320 


254 ± 038 ± 004 


0.57 


0.023 


0.570 


1 576 + 248 + 007 


0.43 


0.057 


0.570 


1 ynq + n osf) + 004 


0.42 


0.100 


0.570 


1 379 ± 226 ± 007 


0.43 


0.167 


0.570 


916 + 151 +0 004 


0.55 


0.336 


0.570 


0.261 ± 0.043 ± 0.002 


0.45 


0.023 


0.795 


1.642 ± 0.244 ± 0.047 


0.45 


0.057 


0.795 


1.581 ±0.241 ±0.011 


0.38 


0.100 


0.795 


1.092 ±0.180 ±0.004 


0.45 


0.167 


0.795 


0.811 ±0.134 ±0.003 


0.54 


0.336 


0.795 


0.228 ± 0.041 ± 0.003 



TABLE XII. CCFR E744/E770 forward differential cross-section for z/^iV jjTix+X at S ~ 209.89 GeV. The forward 

cross-section requires > 5 GeV, and the cross-section values should be multiplied by GpM E /-k. The first error given for 
the cross-sections is statistical and the second systematic. Units are in GeV or GeV^, where appropriate, for the averages of 
the kinematic quantities. 





X 


y 


d<y(v^N^fj, IJ.+ X) 1 

dxdy 1^ 


0.85 


0.023 


0.320 


1.546 ±0.234 ±0.005 


0.53 


0.057 


0.320 


1.613 ±0.256 ±0.011 


0.53 


0.100 


0.320 


1.014 ±0.158 ±0.017 


0.50 


0.167 


0.320 


0.619 ± 0.095 ± 0.004 


0.71 


0.336 


0.320 


0.293 ± 0.041 ± 0.003 


0.68 


0.023 


0.570 


1.812 ±0.300 ±0.015 


0.50 


0.057 


0.570 


1.690 ±0.299 ±0.006 


0.49 


0.100 


0.570 


1.529 ±0.271 ±0.035 


0.50 


0.167 


0.570 


0.756 ±0.134 ±0.003 


0.66 


0.336 


0.570 


0.286 ±0.048 ±0.016 


0.72 


0.023 


0.795 


2.422 ± 0.343 ± 0.027 


0.55 


0.057 


0.795 


2.115 ±0.331 ±0.005 


0.55 


0.100 


0.795 


1.689 ± 0.285 ± 0.005 


0.58 


0.167 


0.795 


0.948 ±0.159 ±0.002 


0.76 


0.336 


0.795 


0.328 ± 0.060 ± 0.002 



TABLE XIII. CCFR E744/E770 forward differential cross-section for v^N iTn+X at S ~ 332.70 GeV. The forward 

cross-section requires E^^+ > 5 GeV, and the cross-section values should be multiplied by G%M E/tv. The first error given for 
the cross-sections is statistical and the second systematic. Units are in GeV or GeV^, where appropriate, for the averages of 
the kinematic quantities. 
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v2 

A. 


X 


v 


dxdy 1^ 


0.36 


0.018 


0.355 


300 + 1 38 + 1 04 


0.41 


0.041 


0.355 


470 ± 177 ± 006 


0.37 


0.069 


0.355 




0.41 


0.111 


0.355 


f) ORQ + n 091 + n 088 


0.60 


0.210 


0.355 


0.126 ± 0.042 ± 0.001 


0.24 


0.018 


0.596 


1 1 35 + 378 + 057 


0.29 


0.041 


0.596 


1 229 + 393 + n 039 


0.31 


0.069 


0.596 


569 ± 176 ± 132 


0.32 


0.111 


0.596 


636 + 205 + 015 


0.50 


0.210 


0.596 


0.181 ±0.058 ±0.009 


0.10 


0.018 


0.802 


1.140 ±0.404 ±0.158 


0.14 


0.041 


0.802 


0.875 ± 0.299 ± 0.086 


0.16 


0.069 


0.802 


0.909 ±0.320 ±0.052 


0.30 


0.111 


0.802 


0.872 ± 0.298 ± 0.036 


0.46 


0.210 


0.802 


0.300 ±0.115 ±0.022 



TABLE XIV. CCFR E744/E770 forward differential cross-section for i/^iV /j,' X a,t E 87.40 GeV. The forward 

cross-section requires E^- > 5 GeV, and the cross-section values should be multiplied by G'pME/tt. The first error given for 
the cross-sections is statistical and the second systematic. Units are in GeV or GeV^, where appropriate, for the averages of 
the kinematic quantities. 





X 


y 


d<y(v^N~^IJ,+ H X) 1 

dxdy 1^ 


0.32 


0.018 


0.355 


0.944 ± 0.371 ± 0.043 


0.48 


0.041 


0.355 


1.116 ±0.386 ±0.085 


0.37 


0.069 


0.355 


0.854 ± 0.299 ± 0.018 


0.40 


0.111 


0.355 


0.570 ±0.191 ±0.004 


0.57 


0.210 


0.355 


0.234 ±0.077 ±0.017 


0.23 


0.018 


0.596 


1.389 ±0.489 ±0.094 


0.37 


0.041 


0.596 


1.542 ± 0.542 ± 0.038 


0.35 


0.069 


0.596 


1.064 ±0.374 ±0.014 


0.35 


0.111 


0.596 


0.718 ±0.256 ±0.010 


0.50 


0.210 


0.596 


0.209 ± 0.071 ± 0.005 


0.11 


0.018 


0.802 


1.699 ± 0.536 ± 0.243 


0.15 


0.041 


0.802 


1.728 ±0.592 ±0.134 


0.27 


0.069 


0.802 


1.734 ± 0.580 ± 0.065 


0.26 


0.111 


0.802 


0.832 ± 0.287 ± 0.020 


0.32 


0.210 


0.802 


0.217 ± 0.079 ± 0.009 



TABLE XV. CCFR E744/E770 forward differential cross-section for v^N ^ ij+ iT X aX E 160.52 GeV. The forward 

cross-section requires E^- > 5 GeV, and the cross-section values should be multiphed by G%ME/tt. The first error given for 
the cross-sections is statistical and the second systematic. Units are in GeV or GeV^, where appropriate, for the averages of 
the kinematic quantities. 
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v2 

A 


X 


V 


dCT(i;^JV-*M+M X) 1 
dxdy 


0.67 


0.018 


0.355 


2 221 + n 714 + ni 3 


0.54 


0.041 


0.355 


1 165 ± 394 ± 015 


0.53 


0.069 


0.355 


1 385 + n 451 +0 012 


0.48 


0.111 


0.355 


n 583 + n 1 89 + 005 


0.68 


0.210 


0.355 


0.152 ± 0.049 ± 0.002 


0.45 


0.018 


0.596 


2 774 + 996 + 1 20 


0.44 


0.041 


0.596 


1 837 + 673 + 054 


0.43 


0.069 


0.596 


1.319 ± 0.501 ± 0.020 


0.43 


0.111 


0.596 


1 078 + 415 + 005 


0.58 


0.210 


0.596 


0.234 + 0.087 + 0.006 


0.17 


0.018 


0.802 


2.710 + 0.796 + 0.280 


0.28 


0.041 


0.802 


1.567 + 0.533 + 0.061 


0.39 


0.069 


0.802 


0.891 +0.334 + 0.033 


0.46 


0.111 


0.802 


0.643 + 0.233 + 0.006 


0.67 


0.210 


0.802 


0.276 + 0.103 + 0.021 



TABLE XVI. CCFR E744/E770 forward differential cross-section for v^,N ii+ fi' X at E ^ 265.76 GeV. The forward 
cross-section requires E^^- > 5 GeV, and the cross-section values should be multiplied by G^ME/tt. The first error given for 
the cross-sections is statistical and the second systematic. Units are in GeV or GeV^, where appropriate, for the averages of 
the kinematic quantities. 



Evis 


Nx>0.5 


MC 




Ns 


Ntot 


M 

M6>5 


34.8-128.6 


1 


3.4 


2.1 


1.3 


688 


0.67 


128.6-207.6 


4 


3.4 


1.9 


1.5 


528 


0.75 


207.6-388.0 


2 


3.4 


2.1 


1.3 


238 


0.78 



TABLE XVII. High-x events using the anti-neutrino data sample, 
for xvis > 0.5, AIC is the Monte Carlo prediction for xvis > 0.5, A'', 
Ntot is the total number of dimuon events, and is the Monte Carlo correction in Eq. IBI 



Evis is in GeV, Nx>o.5 is the number of observed events 
./fc and Ns axe the vr/fc decay and s contributions to A'lC, 



Evis 


Nx>0.5 


MC 


N,/k 


N, 


Ntot 


M 

Me>5 


36.1-153.9 


65 


53.39 


11.6 


2.0 


2304 


0.64 


153.9-214.1 


42 


53.44 


14.8 


1.8 


1598 


0.75 


214.1-399.5 


60 


53.38 


17.4 


2.2 


1201 


0.78 



TABLE XVIII. High-x events using the neutrino data sample. Evis is in GeV, Nx>o.5 is the number of observed events for 
Xvis > 0.5, MC is the Monte Carlo prediction for xvis > 0.5, N^/^ and Ns are the n/k decay and s contributions to MC, 



Ntot is the total number of dimuon events, and 



is the Monte Carlo correction in Eq. Bl 
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FIG. 1. Dimuon production in i/-nucleon DIS from scattering off a strange or down quark (LO QCD charm production). 
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FIG. 2. Schematic of the NuTeV beamline. 
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FIG. 3. Typical dimuon event. 
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FIG. 4. Charged-current Monte Carlo distributions (histograms) compared to data (points). Top: energy of the primary 
muon. Middle: energy of the primary muon entering the toroid. Bottom: hadronic energy. All energies are in GeV. Neutrino 
mode is shown on the left, anti- neutrinos on the right side of the plot. The per degree of freedom is shown for each 
distribution. 
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FIG. 5. Charged-current Monte Carlo distributions (histograms) compared to data (points). Top: neutrino energy in GeV 
Evis = Eifi + E2fi + Ehad- Middle: event vertex position in x plane. Bottom: event vertex position in y plane. Neutrinos are 
shown on the left, anti- neutrinos on the right. The per degree of freedom is shown for each distribution. 
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FIG. 6. X, z, and E distributions for dimuons (points) compared to Monte Carlo (histogram). Neutrinos are shown on the 
left, anti-neutrinos on the right. These distributions are used in the log-likelihood fit. The crosses show the Monte Carlo nK 
background component, and the stars the strange sea contribution. 
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FIG. 7. Dimuon data kinematic distributions (points) not used in the fit: in GeV^, i5^2 in GeV, and y, from top to 
bottom. Neutrinos are shown on the left, anti-neutrinos on the right. The total Monte Carlo prediction is represented by the 
histogram, while the crosses show the Trif background component, and the stars the strange sea contribution. 
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0.126 0.2519 0.3779 0,126 0.2519 0.3779 0.126 0.2519 0.3779 

FIG. 8. (J2fi{x) from NuToV neutrinos, for various Ev — y bins in units of chargcd-currcnt a. The cross-section extracted 

using the BGPAR model in the Monte Carlo is shown in squares, the circles correspond to extraction using the CTEQ model, 
and the triangles to the GRV model. The curves show the model prediction for GRV (dashed), CTEQ (dotted), and BGPAR 
(solid) after the models have been fit to the data. 
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0,1 0,2 0,1 0.2 0,1 0,2 

FIG. 9. (72n{x) from NuTeV anti-neutrinos, for various Ei, — y bins in units of charged-current a. The cross-section extracted 

using the BGPAR model in the Monto Carlo is shown in squares, the circles correspond to extraction using the CTEQ model, 
and the triangles to the GRV model. The curves show the model prediction for GRV (dashed), CTEQ (dotted), and BGPAR 
(solid) after the models have been fit to the data. 
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0.0001 0,1321 0.264 0.396 0.0001 0,1321 0,264 0,396 0.0001 0,1321 0.264 0.396 

FIG. 10. 02fi{x) from CCFR neutrinos, for various Ei, — y bins in units of charged-currcnt a. The cross-section extracted 

using the BGPAR model in the Monte Carlo is shown in squares, the circles correspond to extraction using the CTEQ model, 
and the triangles to the GRV model. The curves show the model prediction for GRV (dashed), CTEQ (dotted), and BGPAR 
(solid) after the models have been fit to the data. 
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0.1 0.2 0.1 0.2 0,1 0.2 

FIG. 11. a2n{x) from CCFR anti-neutrinos, for various E^ — y bins in units of charged-current a. The cross-section extracted 

using the BGPAR model in the Monte Carlo is shown in squares, the circles correspond to extraction using the CTEQ model, 
and the triangles to the GRV model. The curves show the model prediction for GRV (dashed), CTEQ (dotted), and BGPAR 
(solid) after the models have been fit to the data. 
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FIG. 12. Effective DOF and x Bts a function of coverage area C. The open symbols (dotted curve) show the x obtained by 
fitting the cross-section table, while the stars (solid curve) show the effective DOF. The open squares on both curves show the 
result obtained by using the same grid on both smeared and unsmeared variables. 
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